RENORMALIZATION OF SELF-CONSISTENT 
^-DERIVABLE APPROXIMATIONS* 



H. VAN HEES 

Fakultdt fur Physik 
Universitdt Bielefeld 
Universitdtsstrafie 25 
D-33615 Bielefeld 
E-mail: hees@physik.uni-bielefeld.de 

J. KNOLL 

Gesellschaft fur Schwerionenforschung 
Planckstrafle 1 
D-64291 Darmstadt 
E-Mail: j.knoll@gsi.de 



Within finite temperature field theory, we show that truncated non-pcrturbativc 
self-consistent Dyson resummation schemes can be renormalized with local vac- 
uum counterterms. For this the theory has to be renormalizable in the usual 
sense and the self-consistent scheme must follow Baym's <E>-derivable concept. Our 
BPHZ-renormalization scheme leads to renormalized self-consistent equations of 
motion. At the same time the corresponding 2PI-generating functional and the 
thermodynamic potential can be renormalized with the same counterterms used 
for the equations of motion. This guarantees the standard ^-derivable properties 
like thermodynamic consistency and exact conservation laws also for the renor- 
malized approximation schemes. We give also a short overview over symmetry 
properties of the various functions defined within the 2PI scheme for the case that 
the underlying classical field theory has a global linearly realized symmetry. 



1. Introduction 

Describing hot and dense systems of strongly interacting particles, one 
is led to the use of dressed propagators within non-perturbative Dyson- 
resummation schemes. Especially this becomes unavoidable if damping 
width effects are significant for the physical situation in question. 



* based on a talk presented at the conference "progress in noncquilibrium greens func- 
tions, dresden, germany, 19.-22. august 2002" 
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Based on functional formulations by Luttinger and Ward 1 and Lee and 
Yang 2 Baym and Kadanoff 3 studied a special class of self-consistent Dyson 
approximations, which later was reformulated in terms of a variational prin- 
ciple, defining the so-called ^-derivable approximations 4 . The variational 
principle, applied to approximations of the ^-functional, leads to closed 
coupled equations of motion for the mean field and the propagator, which 
guarantee the exact conservation of the expectation values of conserved 
currents and thermodynamical consistency, since at the same time the ap- 
proximated ^-functional is an approximation of the thermodynamic poten- 
tial. 

Later this concept was generalized to the relativistic case and rederived 
within the path integral formalism by Cornwall, Jackiw, and Tomboulis 5 . It 
is no formal problem to extend this formulation to the general Schwingcr- 
Keldysh real-time contour 6 ' 7 and thus to generalize the concept to non- 
equilibrium problems. 

Here we discuss the problem, how to renormalize the equations of mo- 
tion, derived from ^-derivable approximations for relativistic quantum field 
theories. Generalizing the work of Bielajew and Serot 8,9 we show that any 
<f>-derivable approximation of a perturbatively renormalizable theory is also 
renormalizable in the usual sense. Further we prove that at finite tempera- 
ture only temperature-independent counterterms are necessary to give finite 
equations of motion 10 . The counterterms can be interpreted as renormal- 
ization of the wave function and the vacuum parameters of the quantum 
field theory, like the particle masses and the coupling constants. 

Further we demonstrate the possibility to treat numerically ^-derivable 
approximations with generic two-point contributions to the self-energy be- 
yond pure gap-equation approximations with "tadpole self-energies" ) giving 
rise to a finite in-medium damping width of the involved particles 11 . 

Another important question is whether the approximations respect un- 
derlying symmetries of the classical action functional. Contrary to pertur- 
bation theory in general the solution of the ^-derivable equations of motion 
violates the Ward-Takahashi identities of symmetries. This was discussed 
first by Baym and Grinstein 12 on the example of the 0(iV)-symmetric linear 
cr-model. The reason can be traced back to a violation of crossing symme- 
try for approximations of the ^-functional: The solution of the ^-derivable 
equations of motion is equivalent to a resummation of the self-energy in any 
order of the expansion parameter of the ^-functional approximation. This 
involves intrinsically the resummation of higher vertex functions, which is 
incomplete, because certain channels are missing, being taken into account 
only by approximations of the $ functional at higher orders. 
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We show that this can partially be cured by defining a non-perturbative 
approximation to the usual effective quantum action. This admits one 
to define vertex-functions which fulfill the Ward-Takahashi identities of 
the symmetry. These crossing symmetric vertex-functions are defined by 
equations of motion which solution is equivalent to a further resumma- 
tion of the channels, missing intrinsically in the ^-derivable self-energy 
resummation 13 . 



2. The 2PI generating functional 

We start with the defining path integral for the two-particle (2PI) irre- 
ducible quantum action for the state of thermal equilibrium. For this case 
we use the Schwinger-Keldysh closed real-time path extended by an imag- 
inary part making use of the fact that the unnormalizcd thermal density 
operator exp(— /3H), with (3 denoting the inverse temperature of the system 
and H its Hamilton operator, can be included within the path integral as a 
time evolution parallel to the imaginary axis (see Fig. 1). We consider the 
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Figure 1. The Schwinger-Keldysh closed time path modified for the application to ther- 
mal equilibrium quantum field theory. 



local relativistic renormalizable quantum field theory for one scalar field 
with the dynamics defined by the classical action 



S[<f>] 



d(l) 



c 



(1) 



Here and in the following J c d(123 . . .)/i23... denotes an integral over a func- 
tion / of space-time arguments xi,X2, ■ ■ ■■ The time variable is assumed to 
be defined along the contour depicted in Fig. 1. 
The generating functional is given by 

rb,G]=5M + ^Tr(M 2 G- 1 ) + ^ / d(12)Z? 12 1 (Gi 2 - D12) + &[<p, G] (2) 
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with 

°a = P^- P) 

12 8ipi5ip 2 

In terms of diagrams the functional <& consists of closed two-loop dia- 
grams, which are built with lines representing exact propagators G and 
point-vertices with respect to the field <f> derived from the action S[<p + <f>]. 

The equations of motion are determined by the stationary point of the 
functional (2): 

— 1 - (4) 

Sipi SG12 

Using (2) these equations of motion read 



1 — ~ ~n d ^ 2 >^ Gl ' 2 ' _ 1 — ' 



(5) 



Sia := D£ - = 2i^-. (6) 

From the latter equation we see that the derivative of $ with respect to G 
gives the exact self-energy of the theory at presence of a mean field ip, which 
in turn is determined from (5) self-consistently. Since lines in diagrams 
contributing to an expansion of $ stand for exact Green's functions G all 
these lines must not contain any self-energy insertions, i.e., the self-energy 
is represented as the sum of all skeleton diagrams. Since the derivative of a 
diagram with respect to G means to open any line contained in it and then 
adding all the so obtained diagrams, the functional $ consists of all closed 
two-particle irreducible diagrams with at least two loops. 

A ^-derivable approximation is defined as the truncation of the func- 
tional $ to a finite (e.g., coupling-constant or ^-expansion) or an explicitly 
resummable infinite subset of 2PI diagrams. The mean field ip and Green's 
function G are then determined by the self-consistent closed equations of 
motion (5-6). 



3. Renormalization of </> 4 -theory 

For the renormalization of self-consistent approximation schemes we use 
the Bogoliubov-Parasiuk-Hepp-Zimmermann (BPHZ) renormalization de- 
scription. As an example we treat </> 4 -theory in the phase of unbroken 
Z2-symmetry, i.e., we set m 2 > 0. Then as well at zero as at finite tem- 
perature we have ip = as the unique solution of the equations of motion. 
Only the self-consistent propagator has to be determined. 
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3.1. Renormalization at T — 

At T = the only difference to the perturbative treatment of the renormal- 
ization problem is that we have to apply it to diagrams with lines standing 
for self-consistent propagators instead for free ones. It is also clear that 

we can restrict ourselves to the { }-part (i.e., the time ordered part) of 

the real-time contour since in the vacuum case the time-ordered Green's 
function is identical to the retarded (advanced) one for positive (negative) 
p°-components. 

The BPHZ renormalization description rests solely on Weinberg's 
power-counting theorem which is independent of the the special form of 
propagators. Thus we only have to show that the self-consistent propaga- 
tors of <i>-dcrivable approximations belong to the class of functions with the 
asymptotic behavior 0[(Z 2 ) _1 (ln^ 2 )^] for large momenta I 2 , where (3 is a 
constant. Assuming that this is the case Weinberg's theorem tells us that 
a connected truncated diagram 7 with E external lines has a superficial 
degree of divergence S(-f) = 4 — E. So due to the <I>-derivable equations 
of motion the self-consistent self-energy shows an asymptotic behavior like 
0[p 2 (lnp 2 ) /3 ]. So starting an iteration for the self-consistent propagator 
with the perturbative propagator (and provided this iteration converges) 
we can conclude that indeed the propagator is of the usual asymptotic 
behavior. 

The BPHZ renormalization technique aims at the construction of the 
integrand of the renormalized integral without using an intermediate step 
of regularization. If a diagram is divergent without proper divergent sub- 
diagrams it is sufficient to subtract the Taylor expansion of the integrand 
with respect to the external momenta up to the order given by the super- 
ficial degree of divergence which is in our case 4 — E, with E denoting the 
number of external legs. 

For renormalization theory it is crucial that the same holds true for dia- 
grams which contain divergences from proper subdiagrams, if the according 
subdivergences are subtracted first, even if it contains overlapping diver- 
gences and thus that one needs only local counterterms to the quantum 
action which have the same form as that of the classical action but with 
the infinities lumped into the "bare parameters" rather than the physical 
ones 14 ' 15 . 

The described BPHZ-schcmc chooses the renormalization point for di- 
vergent diagrams at external momenta set to 0. It is clear that by another 
finite renormalization of the same diagrams we can switch to any renor- 
malization scheme appropriate for the application under consideration. In 
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our case of (/> 4 -theory we choose the on-shell renormalization scheme, which 
defines the mass parameter m to denote the physical mass of the particles. 
We use the following on-shell renormalization conditions: 

£ (vac V = m 2 ) = 0, d p2 £( vac V = ™ 2 ) = 0, (7) 
r^( S ,t,u = 0) = ± (8) 

Here s, t, u are the usual Mandelstam variables for two-particle scattering, 
p is the external momentum of the self-energy and m 2 is the (renormalized) 
mass of the particles due to the renormalization conditions (7). The second 
condition defines the wave-function normalization such that the residuum 
of the propagator at p 2 — m 2 is equal to unity. 



3.2. Numerical calculation 

To illustrate the abstract considerations of the previous section we show how 
to solve the self-consistent ^-derivable Dyson equation for the 4 model in 
the case of unbroken Z2-symmetry, i.e., for ip = 0. We take into account 
the ^-functional up to three-loop order: 

$ =CO+^ (9) 

For the self-energy we find from Eq. (6) 



+ -< >- (io) 



The main numerical problem is that it is of course not possible to integrate 
directly the renormalized integrands of the self-energy diagrams depicted 
in Eq. (10) because of the on-shell poles of the propagator. Instead we use 
its Lehmann spectral representation 



G< vac V) = / 
Jo 



d(m 2 )ImG( vac H™ 2 ) 



p 2 —\rj 



(11) 



where r] denotes a small positive number to be taken to + in the sense of 
a weak limit after performing the loop integrals. 
First we calculate the one-loop function 

£ (rcs V) = i >0< = 1 j (0f G ^ C)[{l + l) 2 ]G (vac) (l 2 ) (12) 
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which appears as a subdiagram contained in the "sunset diagram" in Eq. 
(10). Here and in the following we use dimensional regularization to give 
the un-renormalized integrals a definite meaning. At the end of the BPHZ 
subtraction procedure we can let d = 4. The loop function (12) is logarith- 
mically divergent and has to be subtracted such that L^ cn \q 2 = 0) = 
due to the renormalization conditions (7). 

Using the spectral representation (11) the renormalized loop function 
can be expressed with help of a kernel K[ rcn ^ : 



Jo ^ Jo n 

x Im G (vac) (ml) Im G (vac) (mp 



(13) 



"2J- 

The renormalized kernel can be calculated analytically with help of stan- 
dard formulae of perturbation theory (for details see 11 ). 

Due to the renormalization conditions the tadpole contribution to the 
self-energy is canceled. For the remaining sunset-diagram we can use the 
dispersion relation for £( ron ) and G^ vac ^ to define a kernel K2 such that the 
renormalized self-energy reads 

r°° Am 2 r°° r\m 2 

S (v a c V)= / dm,/ ™*K™tf,n%,mD 

Jim* K J Q IT ~ (14) 

xImL (rcn) (m^ImG (vac) (m^). 

For the numerical calculation one has to take into account that ImG' vac ' 
contains the pole contribution oc S(p 2 — m 2 ) which has to be treated ex- 
plicitly in both formulae (13) and (14). The remaining integrals over the 
m\ are relatively smooth finite integrals which can be done with help of 
a simple adaptive integration algorithm. We used an adaptive Simpson 
algorithm to solve the equations (13) and (14) iteratively. 

As turns out for the vacuum case the main contribution comes from 
the pole terms such that even for high coupling constants the perturbative 
and self-consistent result lie on top of each other (see Fig. 2). The reason 
is that in our on-shell scheme the threshold for the imaginary part of the 
self-energy is at 9m 2 . 



3.3. Renormalization at T > 

Now we show that the renormalization at T > can be done with the same 
temperature independent vacuum counterterms as were necessary to render 
the vacuum proper vertex functions finite. Thus in complete analogy to 
the well-known result of perturbative finite-temperature renormalization 
theory the renormalized theory is completely defined at T — 0. There 
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Figure 2. 
bative and 
at s = 9m? 



Real (left) and imaginary part (right) of the sunset self- energy. The pertur- 
the self-consistent result lie on top of each other due to the large threshold 



is no ambiguity for "in-medium modifications" of coupling constants from 
renormalization. 

We expand the finite-temperature self-energy around the self-consistent 
solution of the same <I>-derivable approximation at T = 0: 



v — v( vac ) 

^12 — ^12 



y(0) 
^12 



iW 



(15) 



Here 12^°' is the renormalized vacuum self-energy calculated in the previ- 
ous section. The second and third terms in Eq. (15) contain the in-matter 
parts of the self-energy. Thereby is the part of the self-energy which 
arises as the linear part from a functional power expansion with respect to 
the Green's function around the vacuum Green's function: 



-iE 



(0) 
12 



d(l'2') 



5G V 



p(mat) 
U 1'2' 



T=0 



. ir (4) 



(16) 



The wavy line stands for the "matter part" of the Green's function G , ^ mat - ) = 
G-G< vac ). As we shall see below herein we have to understand only the 
diagonal part of the vacuum propagator within the momentum-space matrix 
formalism. The four-point kernel I^ 4 ) is a four-point function represented 
by a particular set of self-energy subdiagrams consisting of pure vacuum 
lines, defined by 



-ir 



(4) 



SE 



12 



6G V 



= -2i 



S 2 $ 



(17) 



T=0 



Its "diagonal part", i.e., with all time arguments placed on one side of the 
real-time contour, defines a vacuum renormalization part which is of super- 
ficial degree of divergence 0. Thus we can conclude that the diagonal part of 
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Q( mat ) ^ f momentum power —4, so that closing I^ 4 ) with a wavy G^ mat - ) - 
line yields another logarithmic divergence, even when the pure vacuum part 
I^ 4 ) is renormalizcd. On the other hand the off-diagonal parts of G( mat ) 
contain ^-functions and Bose-Einstein-distribution factors which lead to 
convergent temperature dependent integrals which we are not allowed to 
subtract. Thus from Weinberg's power-counting theorem we can conclude 
that the divergent part of S^°\ in the following called £(°' dlv ), accounts 
for all terms of momentum power and consequently E^ r ' is of divergence 
degree —2 and thus finite after subtracting vacuum subdivergences. 

So we are left with the task to renormalize the last loop integral from 
closing the r( 4 )-diagram with a G( mat '-line. For this purpose due to our 
discussion above we have to split the full propagator as follows 



1(jt12 — l(jr 12 



i J d(l'2')G^ c) S^ iv) G^ c) 



(r) 



_;p(4,vac) 



iG 



(18) 



Using Eq. (18) shows that £( - dlv ), represented by the second diagram in 
(18), fulfills the equation of motion 



^(Cdiv) 
J 12 



d(i'2')r&£$ I / d(i''2'0GiT^si"; 2 T j G^ a 2 7 + G^, I (19) 



which is linear in G^ r '. Thus Eq. (19) is solved by the ansatz 



-iS 



(0,div) 
12 



— iA( vac ) 



(20) 



1 2 

where the vacuum four-point function A^ vac ^ fulfills the Bethe-Salpeter equa- 
tion 

.(vac) _ F (4,vac) . / , (% . - f ,xp(4,vac) ^(vac) ^(vac) . n 

A l2,l'2' - 1 12, VI' + 1 J c <H«J45bjl 12 34 G 35 G 46 A 56) i/ 2 '- (21) 

Once this logarithmically divergent vacuum subdivergence is renormalized 
also £(°> dlv ) is finite since is falling off with momentum power —6. 

For the renormalization of the four-point function we note that the 
momentum-space version of (21) reads 



A (vac) ^ g j = r (4,vac) ^ q) + i J -jL^T^'^ (p, I) [G^ (I)] 2 (I, (?) 



r (4,vac) fa q )+ij di A (vac) fa Q [ G (vac) (J)]3 r (4,«) ^ ? J 



(22) 



10 



To renormalize this equation, a detailed BPHZ-analysis uses the fact that 
due to the 2PI-feature of the <& functional r^' 34^ is 2PI relative to any cut 
which separates the space-time point pairs (12) and (34). Thus there is no 
"BPHZ-box" cutting through this Bethe-Salpeter kernel. Thus we can do 
the subtractions at the upper and the lower end of any subdivergence with 
the renormalized BS-kernel leading to the renormalized BS-equation 

A {lca) (p,q) = T^ ren \p,q) 

+i / [r( 4 < rcn ) (p, i) - r( 4 < vac ) (o, 0] [G {vac) (i)} 2 A (icn) (23) 
+i/ I g T A( ron )(o,0[G( vac )(0] 2 [r (4 ' ron) (^'z) - r( 4 - ren )(z,o)]. 

The renormalization of r( 4 ' vac ) itself is straight forward, since it is given 
by a finite set of vacuum diagrams which can be renormalized by the same 
BPHZ-scheme as the perturbative ones. 

For the practical calculation of the self-energy part S^ -* we need only 
the A( rcn )(0,g). Indeed using Eqs. (18) and (20) we find 

£<°>(p) - £(°>(p) - £(°)(p) + S(°)(p) 

= / -0^ [r( 4 - ren ) (p, i) - r( 4 < rcn ) (o, z)]G( mattor ) (/) (24) 

+ /|| I AM ( o i ;)GM(i) 

An example solution for the equations and the comparison with the per- 
turbative approximation is shown in Fig. 3. 

While at T = both, the perturbative and the self-consistent solution, 
show the three-particle threshold at y/s — 3m at finite temperature the 
spectral width smoothes out all threshold structures in the self-consistent 
solution. The growing high-energy tail is related to the decay of a virtual 
particle to three particles. At finite temperature as an additional effect a 
low-energy plateau in Im T, R emerges from in-medium scattering processes 
of real particles. 

The comparison of the self-consistent solutions with the perturbative 
approximation shows counterbalancing effects of self-consistency: The fi- 
nite spectral width, contained in the self-consistent propagator leads to a 
further broadening of the width and a smoothing of the structure as a func- 
tion of energy This is counterbalanced by the behavior of the real part of 
the self-energy which essentially shifts the in- medium mass upwards. This 
reduces the available phase space for real processes. With increasing cou- 
pling strength A a nearly linear behavior of Im E fl with po results, implying 
a nearly constant damping width, given by — ImE^/po- 

While the tadpole contribution always shifts the mass upwards, at 
higher couplings and temperature also the real part of the sunset diagram 
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Pert. Re !2 for 1=30, T=250 MeV Perl. Im £ for 1=30, T=250 MeV 
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Figure 3. Real (left) and imaginary part (right) of the perturbativc (top) and the self- 
consistent self-energy for A = 30, m = 140MeV and T = 250McV. Note that the 
self-energies are multiplied with an factor 100 in these plots! 



becomes significant which can lead to a net down-shift of the mass again. 
This can be seen for the parameter set used for the Fig. 3. 



4. Global symmetries 

If the classical action underlying a quantum field theory is invariant under 
the linear operation of a Lie group G, where the group elements are inde- 
pendent of the space-time argument one speaks of a global symmetry of the 
classical action. Then it can be shown that for each linearly independent 
generator t a e CG there exists a conserved charge Q a which in turn builds a 
basis of the Lie algebra CG as a subalgebra of the canonical Poisson algebra 
of the fields. 

It is a well-known theorem in perturbativc quantum field theory that 
the same holds true for the quantized theory i.e., the effective action is 
symmetric under the same symmetry group as the classical action provided 
that no quantum anomaly destroys the symmetry or a part of it. The n- 
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point Green's functions fulfill certain constraints due to the symmetry which 
we shall call the Ward-Takahashi identities (WTIs) of the symmetry 

For <I>-derivable approximations in general the WTIs of the symmetry 
are violated for the two-point and higher vertex functions although the 
expectation values of the corresponding Noether currents are exactly con- 
served. This can be traced back to a violation of crossing symmetry within 
the self-consistent propagator due to partial resummation of the self-energy 
insertions by means of the Dyson equation of motion. We shall further show 
that the WTIs are fulfilled for vertex functions, which we shall call external. 
These functions arc defined from a non-perturbative effective action which 
is uniquely determined by the underlying ^-derivable approximation. 

We study the symmetry properties of the generating functional (2) of 
a scalar 0(iV)-symmetric quantum field theory. An infinitesimal O(N)- 
transformation reads 8(jy{ = \5x\ {j a y y<\>\ ■ Making use of the assumed 
invariance of the classical action under 0(iV)-transformations one finds by 
standard path-integral analysis 13 that T[ip, G] is an 0(7V)-scalar functional 
when <f> is transformed as a vector and G as a second-rank tensor: 

/ d(l)-£(r»)V? + / d(12)-^[(r a )VG{ 2 fe + (r a )\,G^'] = 0. (25) 
Jc oq>\ Jc 5G{ 2 

Now only for the exact functional the self-consistent Green's function G 
is identical with the exact one. Thus in general only for the exact case 
the self-energy and higher vertex functions fulfill all WTIs of the usual 1PI 
action. 

Generally for a <I>-derivable approximation this equality of vertex func- 
tions does not hold true any longer. For such approximations in general the 
WTIs are violated in loop orders higher than that taken into account for 
the approximation of the functional. The reason is that the solution of the 
equations of motion are equivalent to a certain partial resummation of the 
perturbation series which is not crossing symmetric in the inner structure 
of the diagrams. 

To recover the crossing symmetry we define an effective action functional 
from our approximated 2PI-action as 



r[<p]=T[p,G[<p]] with ST ^ G ^ 



5G 



= 0, (26) 

G=G[<p] 



i.e., for an arbitrarily given mean field <p we define the propagator G[<p] 
as the solution of the Dyson equation of motion as it is defined by the 
2PI- functional. The approximation to the 1PI functional T[ip] is given by 
insertion of this propagator in the 2PI functional. If one would not ap- 
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ply any approximations f is identical with the usual 1PI quantum action 
functional generating 1PI truncated proper vertex functions. 

The stationary point of this action functional defines the mean field and 
propagator of the ^-derivable approximation since 

st _ ( st[ v ,o\ , / A „ m sr[v,a\ sg{%\ \ 

W-\ *> + L m) % ) a __ m <27) 

From (26) this yields 



ST = f ST[<p,G\ \ 

<V V 8<P )g=&[ v ] 



(28) 



Thus the stationary point (p of T is identical to the mean field of the <f>- 
derivable approximation and together with (26) this means that G[ip] is the 
solution of the Dyson equation from the same ^-functional. 

For a ^-derivable approximation T defines a non-perturbative approxi- 
mation to this functional and can be used to derive approximations for the 
proper vertex functions: 

bip\6ipl ■ ■ ■ 

Especially is an approximation for the inverse propagator, which we 
call the external propagator to be distinguished from the self-consistent 
propagator. This external propagator fulfills the usual WTI but is not 
identical with the self-consistent or internal propagator. 

The vertex functions (29) are by definition crossing symmetric. Now the 
symmetry property (25) by construction holds also true for the approxi- 
mated functional and thus from (26) we see that T is an 0(7V)-scalar 
functional: 

d(l)^(r%Vi'=0. (30) 

This symmetry property contains all WTIs for the external vertex functions. 
Especially for the external propagator, defined by 



l G extA.J,2fc 



By taking the functional derivative of (30) we obtain 



(31) 



|d(l)(C)i^(r7/i=0. (32) 
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For a translationally invariant state the Fourier transform of (32) with 
respect to (x\ — x 2 ) reads 



0. 



(33) 



It is clear that in this situation ip is a constant due to translation invariance. 
If it is not the symmetry is spontaneously broken. Since (M 2 )jk is the 
mass matrix of the particles described as the excitations of the field around 
the mean field (p Eq. (33) tells us that the N — 1 fields perpendicular to 
the direction given by the solution (p are massless, the Nambu-Goldstone 
bosons of the symmetry. 

To calculate the external propagator explicitly we apply (26) in (31) to 
obtain 



( G oxt)lj,2*; 



Taking the derivative of the identity 



c V<5G^: Srf 



- <p=(f, G—G[(p\ 



d(2')(G- 1 )i J , 2 ' fc 'G^' 2 fc = S^S* 



with respect to the field we get 
'8G- 1 



d(l'2') 



5<p 3 



G~ k k i ( ri— 1\ 
2'2 + \G )lj,2'k' 



5Gt 



k'k 



2'2 



With help of (6) the three-point function 



,(») 



-y,2fe; 3 ; - ^ 
can be expressed as the solution of the BS-equation 



(34) 



(35) 



(36) 



(37) 



^lj,2k;3l — ^lj,2k;3l~^ d ( 3 ' 4 ' 3/ '4'0 r i^2fc;3'i'4'm'G34''G4™4™ A^„ i4 » ro // ;3 ; , (38) 

which in terms of diagrams can be depicted as 
iA< 3 ) 



while its kernels are determined by the ^-functional: 




r (3) 
l?,2fc;3Z 



r (4) 

L lj,2k;3lAm 



S 3 S[ V ] 2 . s 2 Hp,g] 



6G{ k 2 6G l ™ 



G=G[<p], p—(p 



(40) 



G—G[if], (p—p 
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With these definitions using (36 and (37) the external propagator (34) can 
be written as 



(£*ext)lj,2fc — 



2 Jc Sip^dip^Sifl Sip*' <Vi<fy'2 

i jf d(3'4'3"4")rg, 4 , fe , ;1 ..Gj'^G^As^,, 4 „ fc „ 



2A- 



(41) 



In graphical terms this equation looks as follows 



(42) 



It is clear that both, the BS-equation (38) and Eq. (41), have to be rcnor- 
malized. This is done again with help of the above explained BPHZ tech- 
niques. Again all counterterms turn out to be independent of the tem- 
perature and consistent with those needed to renormalize the underlying 
self-consistent equations of motion (for more details of this renormalization 
see 13 ). 

Looking at the diagrams (39) and (42) it turns out that the BS-equation 
(38) provides exactly the resummation of the channels missing inside the 
self-consistent approximation to the self-energy. This restores both the 
intrinsic crossing symmetry and the underlying 0(A r )-WTI for the external 
self-energy. 

As an example we show results for the lowest order approximation which 
is the Hartree approximation for the self-consistent self-energy, leading to 
a constant effective mass shown in Fig. 4. Clearly Goldstone's theorem 
is violated, since the pion mass is different from in both broken phases. 
In this case the external self-energy is obtained by a bubble resummation, 
corresponding to a Random phase approximation (RPA). Since at the same 
time it is the second functional derivative of the effective lPI-action func- 
tional it provides a stability criterion for the Hartree solution, which is of 
course only a stable solution, i.e., a minimum of the effective potential, if 
the mass matrix is positive semidcfinite. It turns out that the solution, 
denoted as "broken phase 2" in Fig. 4 is unstable. In Fig. 5 the effective 
mass obtained from the external propagator, is shown. Indeed Goldstone's 
theorem is fulfilled, as to be expected from our analysis. From another 
point of view this result was also obtained by Aouissat and Belacem 16 . 
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Self-consistent a-mass at finite temperature 
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Self-consistent it-mass at finite temperature 
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Figure 4. The solutions for the Hartree approximation 
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Figure 5. The effective external masses at a temperature of 150MeV. The effective ex- 
ternal 7r-mass indeed vanishes at po = p = as predicted from Goldstone's theorem. The 
spectral function of the tr-meson shows that at high temperatures its strength becomes 
more peaked and the maximum shifted to lower momenta than at T = 0. 



5. Conclusions and outlook 

We have shown that any self-consistent Dyson resummation can be renor- 
malized with counter terms that are independent of temperature provided 
it is realized as a ^-derivable approximation. The proof is based on the 
applicability of Weinberg's power counting theorem and the BPHZ sub- 
traction scheme. These techniques provide the possibility to extract both, 
the explicite and hidden divergent vacuum subdiagrams, and to subtract 
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the divergences leading to a coupled set of finite equations of motion for 
the mean fields and the Green's functions. 

In the second part a detailed analysis of symmetry properties has been 
given. In general the ^-derivable approximation schemes violate Ward- 
Takahashi identities for two-point and higher Green's functions since in- 
trinsically crossing symmetry is violated at orders of the expansion param- 
eter higher than that taken into account for the 2PI functional. Also the 
recently observed deviation from the correct renormalization group flow of 
the coupling constant can be traced back to this incomplete resummation 
provided by the self-consistent equations of motion 17 . 

It was further shown that the symmetries arc recovered by defining 
a non-perturbative approximation to the effective quantum action by the 
stationary point of the 2PI action functional with respect to the Green's 
function at given mean fields. The calculation of the self-energy defined 
with help of this 1PI functional needs the solution of a Bethe-Salpeter 
equation keeping track of the channels that are missing intrinsically in the 
(f>-derivable self-consistent resummation for the self-consistent propagator. 
In this way not only crossing symmetry but also the Ward-Takahashi iden- 
tities for the proper vertex functions related to linearly realized global sym- 
metries are recovered. 

The problems related with the missing crossing symmetry and symmetry 
violations become more serious in the case of local gauge symmetries. In 
general then the <I>-derivable approximations violate important features like 
unitarity, causality and positive definiteness of the probability measure. 
Here only partial solutions of these problems are known, for instance making 
use of hard thermal loop expansion schemes 18 ' 19 . A systematic analysis of 
the violation of gauge invariance has been given recently by Arrizabalaga 
and Smit 20 . It was shown by Dcnncr and Dittmaier 21 that there exist non- 
symmetry breaking Dyson resummation schemes within the background 
field gauge formulation. At the time an investigation whether or not this 
scheme is applicable also in the context of the ^-derivable formalism is 
under way. Nevertheless the definition of the self-consistent propagator is 
in this case questionable because of the artificial excitation of unphysical 
gauge field degrees of freedom by the violation of crossing symmetry inside 
the self-consistent diagrams. A first way out of this problem by a projection 
formalism is given by us in 22 . 

For the study of non-equilibrium situation the ^-functional becomes 
the only systematic approximation scheme which obeys conservation laws. 
Recently numerical studies of the equations of motion for out of equilibrium 
were undertaken (see 23 > 24 > 25 > 26 and citations therein). 
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The ^-derivable approximations are also a starting point for the deriva- 
tion of transport equations from quantum field theory, ensuring consistency 
conditions like detailed ballance, Boltzmann's H-theorem and conservation 
laws. Especially it becomes important for finding consistent transport equa- 
tions for particles or resonances with a finite mass width 27 ' 28 ' 29 ' 30 . 
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